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S-nitrosoglutathione reductase (GSNOR) is believed to modulate effects of reactive oxygen 
and nitrogen species through catabolism of S-nitrosoglutathione (GSNO). We combined 
bioinformatics of plant GSNOR genes, localization of GSNOR in Arabidopsis thaliana, 
and microarray analysis of a GSNOR null mutant to gain insights into the function and 
regulation of this critical enzyme in nitric oxide (NO) homeostasis. GSNOR-encoding genes 
are known to have high homology across diverse eukaryotic taxa, but contributions of 
specific conserved residues have not been assessed. With bioinformatics and structural 
modeling, we show that plant GSNORs likely localize to the cytosol, contain conserved, 
solvent-accessible cysteines, and tend to be encoded by a single gene. Arabidopsis 
thaliana homozygous for GSNOR loss-of-function alleles exhibited defects in stem and 
trichome branching, and complementation with Green fluorescent protein (GFP) -tagged 
GSNOR under control of the native promoter quantitatively rescued these phenotypes. 
GSNOR-GFP showed fluorescence throughout Arabidopsis seedlings, consistent with 
ubiquitous expression of the protein, but with especially high fluorescence in the root 
tip, apical meristem, and flowers. At the cellular level we observed cytosolic and nuclear 
fluorescence, with exclusion from the nucleolus. Microarray analysis identified 99 up- and 
170 down-regulated genes (>2-fold; p< 0.01) in a GSNOR null mutant compared to wild 
type. Six members of the plant specific, ROXY glutaredoxins and three BHLH transcription 
factors involved in iron homeostasis were strongly upregulated, supporting a role for 
GSNOR in redox and iron metabolism. One third of downregulated genes are linked 
to pathogen resistance, providing further basis for the reported pathogen sensitivity 
of GSNOR null mutants. Together, these findings indicate GSNOR regulates multiple 
developmental and metabolic programs in plants and offer insight into putative routes 
of post-translational GSNOR regulation. 
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INTRODUCTION 

In plants, biological processes ranging from leaf stomatal clo- 
sure to auxin perception in roots and pathogen infection involve 
nitric oxide (NO) (NeiU et al., 2002; Floryszak-Wieczorek et al., 
2007; Lozano-Juste and Leon, 2011; Terrile et al, 2012). While 
NO itself is ostensibly active, it is also thought to be transmitted 
to distal targets via low molecular weight S-nitrosothiols (SNOs), 
of which the glutathione (GSH) adduct S-nitrosoglutathione 
(GSNO) is the most abundant (Broniowska et al, 2013; Corpas 
et al., 2013). GSNO can profoundly affect protein activity through 
glutathionylation and nitrosation of cysteines (Romeo et al., 
2002; Giustarini et al., 2005; Zaffagnini et al., 2013), imply- 
ing that cells require mechanisms to spatiotemporally control 
GSNO levels. Catabolism of GSNO by S-nitrosoglutathione 
reductase (GSNOR) is common to eukaryotes and many bac- 
teria and is believed to be responsible for this regulation of 



GSNO levels (Liu et al, 2001; Staab et al, 2008). GSNOR exhibits 
NAD/H-dependent oxidoreductase activity toward a broad spec- 
trum of aliphatic compounds, but its preferred substrates are 
GSNO and S-hydroxymethylglutathione (HMGSH), an interme- 
diate in formaldehyde metabolism (lensen et al., 1998; Achkor 
et al, 2003; Kubienova et al, 2013). While GSNO catabolism has 
been observed with Cu-Zn superoxide dismutase, GSH peroxi- 
dase, xanthine oxidase, and human carbonyl reductase 1 (CRl), 
the former three enzymes merely regenerate NO [summarized 
in Broniowska et al. (2013)], and residues critical to interac- 
tion between CRl and GSH adducts are not conserved in plants 
(Bateman et al, 2008). GSNOR is therefore considered the pri- 
mary catalyst for GSNO catabolism in plants. 

The importance of GSNOR to plant growth, develop- 
ment and stress responses has been highlighted by several 
studies. Lowered GSNOR expression in Arabidopsis thaliana 
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(Arabidopsis), resulting from a null mutation {atgsnorl-3/hot5-2) 
or RNAi, was correlated with higher SNO content and differen- 
tial susceptibilities to pathogens (Feechan et al., 2005; Rusterucci 
et al., 2007). The effect of absence of GSNOR was extended by 
Lee et al. (2008 ) who described a thermotolerance defect that was 
rescued with NO scavengers. Other phenotypes of plants with 
GSNOR mutations include diminished fertility and resistance to 
programmed cell death induced by paraquat, an herbicide that 
elicits robust reactive oxygen species (ROS) production (Lee et al, 
2008; Chen et al, 2009). These concomitant gains and losses 
of function are analogous to consequences of GSNOR inhibi- 
tion in mammals, for which both enhanced carcinogenesis and 
abated severity of inflammatory diseases are observed (Wei et al., 
2010; Sun et al., 2011; Tang et al, 2013). Such pleiotropy suggests 
GSNOR participates in both homeostatic maintenance and biotic 
and abiotic stress responses. 

The evolutionary conservation of GSNOR is high (Liu et al, 
2001), and although the consequences of GSNOR depletion have 
been described at the organismal level for Arabidopsis (Lee et al, 
2008; Chen et al, 2009; Kwon et al, 2012), a molecular etiology 
for the GSNOR loss-of-function phenotype is lacking. Here we 
sought to address the issues of how GSNOR activity could be reg- 
ulated and of what processes are impacted by changes in GSNO 
levels and, therefore, potentially regulated by nitrosation or glu- 
tathionylation of protein effectors. We searched for conserved and 
unique features of plant GSNOR proteins, localized GSNOR at 
the tissue and cellular levels, and measured global changes in 
the transcriptome of an atgsnor/hotS nuU mutant in Arabidopsis. 
Our data demonstrate that most sequenced green plant genomes 
are predicted to encode a single copy of GSNOR characterized 
by a high content of positionally-conserved cysteines. GSNOR is 
found in the cytosol and nucleus throughout the plant, and is 
thus available to modulate GSNO concentration in most if not all 
cells. Moreover, alterations in the transcriptome of Arabidopsis 
homozygous for the atgsnorl-3/hot5-2 nuU allele exhibited dys- 
regulated expression of pathogen response and calcium signaling 
genes, but higher expression of a subset of glutaredoxin (GRX)- 
encoding genes. Together, these data suggest GSNOR facilitates 
multiple homeostatic and stress adaptation processes in green 
plants. 

MATERIALS AND METHODS 

AMINO ACID SEQUENCE ALIGNMENTS AND ANALYSIS 

Predicted genes encoding GSNOR from green plants were 
retrieved from NCBI Genbank and Phytozome v9.1 (Goodstein 
et al, 2012) using Arabidopsis GSNOR (At5g43940, ADH2) 
as the tblastn query. GSNOR copy number was assessed with 
Phytozome and NCBI tblastn algorithms by querying predicted 
GSNOR- encoding genes against genomic reads from a particu- 
lar plant species. Similar hits were considered duplicates if 5' and 
3' intragenic and intron nucleotide sequences were >99% iden- 
tical. Predicted GSNORs were also aligned via ClustalW (Larkin 
et al., 2007) with the Arabidopsis protein most similar to GSNOR 
(alcohol dehydrogenase, AtADHl, Atlg77120), and sequences 
that cladded with AtADHl were culled. Bacterial, metazoan, and 
fungal orthologs discussed in the text were uncovered through 
an NCBI tblastn search with E. coli, Saccharomyces cerevisiae, 



and human GSNOR queries, respectively. ClustalW sequence 
alignments were made with Jalview (Waterhouse et al., 2009). 
Deduced protein sequences were included in Figure 1 if tran- 
script evidence (i.e., RNAseq and/or ESTs spanning the coding 
sequence) was available through Phytozome. Phylogenetic trees 
were drawn from aligned sequences in EvolView (Zhang et al, 
2012). N-terminal targeting peptide searches were performed 
with Predotar (Small et al, 2004) and MITOPROT (Claros, 1995). 
Mitochondrial targeting peptides encoded by GSNOR 5' intra- 
genic regions were identified by six-frame translation of the first 
500 base pairs upstream of the start codon using the intron splice 
rules for Arabidopsis (Hebsgaard et al., 1996). 

GSNOR STRUCTURAL ASSESSMENTS 

Crystal structures of GSNOR from tomato (4DL9), human 
(IMPO), and Arabidopsis (4JJI and 4GL4) were obtained from 
the Protein Data Bank. Graphics were made with PyMOL v 1.6 
(PyMOL). 

PLANT MATERIAL AND GROWTH CONDITIONS 

The Arabidopsis GSNOR nuU mutants, hot5-2 (also known as 
atgsnorl-3) (Col-0 background) and hot5-4 (WS background), 
have been described previously (Lee et al., 2008). Unless otherwise 
indicated, plants were grown in soil in growth chambers on 16h 
days (150 [xMol m"^ s"^ light intensity) and a 21/19°C day/night 
temperature cycle. For analysis of the number of branches pro- 
duced by wild-type and mutant plants, plants were grown as 
above and their height and branch pattern and numbers were 
measured 8 weeks after germination. 

GENERATION AND VISUALIZATION OF GSNOR-GFP FUSIONS IN 
PLANTS 

Green fluorescent protein (GFP) was fused in frame 3' to the 
Arabidopsis GSNOR genomic DNA (including 754 bp and 180 
bp of 5' and 3' UTR, respectively) as follows. Genomic DNA 
was cloned into pENTR/D TOPO (Life Technologies). GFP 
was obtained from pMDC83 (Curtis and Grossniklaus, 2003). 
Hindlll and Xbal sites were added to GFP and GSNOR-encoding 
sequences with primers GGATGCAAGCTTAGTAAAGGAGAAG 
AAC and TCTCTAGATTATTTGTATAGTTCATCCATGC (GFP) 
and primers CGTTGTG TCCTCGATACCAGCAAGCTTGTC 
TCTAGATGACTATATGGGTCCTCTCTGC and GCAGAGA 
GGACCCATATAGTCATCTAGAGACAAGCTTGCTGGTATCG 
AGGAC ACAACG (GSNOR). PGR products were digested with 
Hindlll and Xbal, purified, and ligated. The GSNOR-GFP clone 
was then subjected to PGR with primers GGTACCGAA TTC 
CTAGAGTACAACCTC and TCGAGTGCGGCCGCTAAACT AT 
ATGATTAG to add EcoRl and Notl sites, respectively The PGR 
product was then ligated into _EcoRI/NotI-digested pENTR2B. 
Constructs were cloned into the binary vector pBIB Basta-GWR 
(Gou et al, 2010) with LR Clonase II (Invitrogen), after which 
the sequence was verified. 

The GSNOR-GFP translational fusion was then transformed 
into hot5-2 plants, and homozygous transformants were iden- 
tified by Basta screening and western blots. Three independent 
transgenic lines were analyzed. Whole seedlings and roots were 
imaged 7 days after germination on minimal nutrient medium 
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FIGURE 1 I Continued 

Conserved features of plant GSNOR proteins. Dark blue, liglit blue, and 
uncolored residues, respectively, refer to 100%, >75%, and <75% sequence 
conservation. Black and red dotted horizontal lines demarcate the catalytic 
and NADH-binding domains, respectively, as reported in crystal structure of 
tomato GSNOR (PDB code 4DL9). Residues coordinating structural and 
catalytic zinc atoms are outlined by solid black and red boxes, respectively. 
Red asterisks denote substrate-binding amino acids according to Kubienova 
et al. (2013). Dotted black boxes highlight flexible regions enclosing the active 
site. Red arrowheads (T) indicate positions of the hot5-l and hot5-3 
missense mutations (Lee et al., 2008). Open (V) and closed (T) black 



with 0.5% sucrose. Whole flowers and reproductive structures 
were analyzed using stage 13 or 14 flowers as indicated in the fig- 
ure legend. Images were obtained using an Olympus Fluoview 
FVIOOO confocal microscope, with the exception of the whole 
seed and isolated stamens, which were imaged using conventional 
fluorescence microscopy with a NIKON Eclipse E800 microscope 
equipped with a SPOT camera (Molecular Diagnostic). 

MICROARRAY ANALYSIS 

Wild-type (Col) and hot5-2 Arabidopsis plants were grown on 
soil in a growth chamber on a 12 h light, 21°C/12h dark, 19°C 
cycle for 25 days after germination. Four biological replicates of 
wild-type and hot5-2 leaves were sampled for RNA extraction 
1 h before the end of the light period. A two-color, dye swap 
hybridization was performed on a long-oligonucleotide array 
chip by the Galbraith lab (University of Arizona) according to 
published methods (Zanetti et al, 2005; Zhang et al, 2008). Data 
were analyzed with Robin Version 0.9.6 BETA (Lohse et al, 2010). 
Differentially-expressed genes were identified by two criteria: (1) 
change in expression greater than 2-fold, and (2) a t-test p-value 
<0.01. Normalization of expression data, analysis using a linear 
model, and Benjamini and Hochberg false discovery rate cor- 
rection for multiple comparisons were performed using LIMMA 
(Smyth and Speed, 2003). 

RESULTS 

CATALYTIC AND ZINC-COORDINATING RESIDUES ARE CONSERVED IN 
GSNOR FROM GREEN PLANTS 

To identify potentially novel shared motifs in plant GSNORs, 
we employed the tblastn algorithms of Phytozome (Goodstein 
et al., 2012) and NCBI Gene Expression Omnibus (Altschul et al, 
1997) to search for GSNOR sequences supported by transcrip- 
tional data. Eighteen unique, type III ADH-encoding cDNAs 
were obtained from a variety of monocots, dicots, mosses, and 
protists. Residues near the dimer interface of tomato GSNOR 
(Kubienova et al., 2013), whose mutation compromises thermo- 
tolerance in dark grown Arabidopsis seedlings (Lee et al, 2008), 
are notably found in all plant sequences, while a glycine that when 
mutated to aspartate diminishes GSNOR activity, but confers 
enhanced paraquat resistance (Chen et al, 2009), is present in all 
but two algae (Figure 1). Substrate- and NADH-enclosing clefts 
(Figure 1, black dotted boxes) are identical or contain conser- 
vative substitutions among moss and algal orthologs. Structural 
and catalytic zinc-coordinating residues (Figure 1, solid black and 
solid red boxes, respectively) and substrate-binding amino acids 



arrowheads designate ex-zinc cysteines found in most and all plant 
sequences, respectively. Z mays: maize. S. lycopersicum: tomato. O. sativa: 
rice. R vulgaris: bean. R patens: Physcomitrella. M. truncatula: Medicago. C. 
reinhardtii: Chlamydomonas. S. moellendorffii: Selaginella. G. max: soybean. 
G. raimondii: cotton. R trichocarpa: poplar. R persica: peach. O. lucimarinus: 
Ostreococcus. S. tuberosum: potato. S. cerevisiae: budding yeast. H. 
sapiens: human. Black asterisks: species with genomes predicted to encode 
additional paralogs, but only transcript-supported sequences are shown. 
Accession numbers can be found in Supp. Table 1 and an alignment of 
additional plant sequences not currently supported by transcript data is 
presented in Supp. Figure 1. 



(red asterisks), as reported by Kubienova et al. (2013), are identi- 
cal in all but two predicted proteins, which contain a single lysine 
to arginine, conservative substitution. 

GSNOR is remarkably cysteine rich, with a mole percent 
cysteine of 3.84 % for the Arabidopsis protein, compared to 
the 1.37 % average for aU proteins in the UniProtKB database 
(2013). Because cysteines can serve as key post-translational reg- 
ulatory sites being modified by nitrosation, glutathionylation, or 
reversible oxidation, we analyzed the conservation of the nine 
extra non-zinc-coordinating cysteine residues (ex- zinc cysteines) 
in Arabidopsis GSNOR. Four are in all the transcript-supported 
plant sequences, two are substituted in one, two are substituted 
in two, and one differs in four organisms (Figure 1, solid and 
open black arrowheads), yielding an overall conservation of 93.8 
% (i.e., on average, each cysteine is present in 17 of the 18 
sequences). If additional plant genes are considered for which 
EST support is lacking, ex- zinc cysteine conservation is stiU 91.0% 
(Supp. Figure 1 — see Supp. Table 1 for accession codes). Thus, the 
position and frequency of ex-zinc cysteines are highly conserved 
in plant GSNORs. We also examined the position of ex-zinc 
cysteines in the Arabidopsis GSNOR structure (PDB 4JJI, via 
(PyMOL) and found that three were solvent accessible (Cys-10, 
Cys-271, and Cys-370, Figures 2A-C), two of which — Cys-10 
and Cys-271 — are positionally conserved even in the human 
sequence. The structures of GSNOR from human (IMPO) and 
tomato (4DL9) showed similar solvent exposure of the homol- 
ogous ex-zinc cysteines, suggesting these three residues may have 
conserved functions in regulating GSNOR activity. 

MOST PLANT GENOMES ENCODE ONE GSNOR PROTEIN PREDICTED 
TO BE FOUND IN THE CYTOSOL 

Encoded by a single gene, Arabidopsis GSNOR consists of nicoti- 
namide cofactor-binding and catalytic domains and has two 
primary enzymatic activities — GSNO reductase and HMGSH 
dehydrogenase (Lee et al., 2008; Crotty, 2009; Kubienova et al., 
2013). Interestingly, although GSNOR is found primarily as a 
single-copy gene in most other examined plant species (12 of 15 
transcript-supported organisms, and 35 of 41 green plants ana- 
lyzed), the dicots Populus trichocarpa (poplar) and Gossypium 
raimondii (diploid cotton) and the moss Physcomitrella patens 
(Physcomitrella) are predicted to have two GSNORs. Although 
not yet supported by transcript data, gene families are also 
predicted in Phaseolus vulgaris (bean). Glycine max (soybean), 
and Malus domestica (apple; Supp. Figure 1). Moreover, in 
phylogenetic trees calculated by the average distance method 
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FIGURE 2 I Of nine positionally-conserved ex-zinc cysteines in GSNOR 
tliree are solvent-accessible. Three orientations (rotation angles indicated) 
of a monomer of tfie Arabidopsis GSNOR dimer are shown with solvent 



accessible surface (PDB 4JJI, 1.8 A res., Rfree = 0.223) in yellow. 
Solvent-accessible ex-zinc cysteines are indicated in blue, and the dimer 
interface in gray. Images were made in PyMOL. 



with EvolView (Zhang et al, 2012), paralogs within the same 
species were more similar to one another than to orthologs 
in other species (data not shown). Thus, not only do most 
plant genomes encode one copy of GSNOR, but duplication of 
GSNOR-encoding genes has occurred recently and sporadically 
among plant species. 

Extended N-termini, which could facilitate organelle targeting, 
are present in predicted GSNOR sequences from Pninus persica 
(peach) and Physcomitrella (Figure 1). Therefore, the Predotar 
(Small et al., 2004) program was used to search for putative 
endomembrane, plastid, and mitochondrial targeting peptides in 
all GSNORs. Of 50 input sequences, only one Physcomitrella par- 
alog was predicted to have a mitochondrial targeting peptide. 
Analysis with MITOPROT (Claros, 1995) achieved similar results. 
Intriguingly, 5' intragenic regions of GSNORs from Arabidopsis, 
rice, and Selaginella moellendorfii included cryptic splice sites that 
could give rise to putative mitochondrial targeting peptides, but 
intragenic regions from tomato, potato, Medicago, and O. luci- 
marinus did not. These results suggest GSNOR is a cytosolic 
enzyme in most plants. 

We also noted predicted N- and C-terminal extensions in non- 
transcript-supported orthologs from two strains of Micromonas 
pusilla (an alga) and in apple that are too long to be signaling pep- 
tides (Supp. Figure 1). BLAST searches indicated M. pusilla exten- 
sions were formylglutathione hydrolase (FGH) domains, which 
catalyze the decomposition of S-formylglutathione (FGSH) to 
GSH and formate. FGHs are predicted in many green plants, 
albeit as separate gene products. FGSH is the product of the 
HMGSH dehydrogenase activity of GSNOR and is a key inter- 
mediate in formaldehyde detoxification (Staab et al, 2008). This 
observation suggests that, at least in the case of M. pusilla strains, 
NO and formaldehyde metabolism are intimately linked. The 
C-terminus of an apple GSNOR comprises a 4,5-DOPA dioxyge- 
nase (DOD), an enzyme in the biosynthesis of phenylpropanoid- 
derived betalain pigments that aid in plant defense (Georgiev 
et al, 2010), but transcriptional support has not been obtained 
for this fusion. 



EXPRESSION AND LOCALIZATION OF ARABIDOPSIS GSNOR 

To determine the major sites of GSNOR function, we created 
a translational fusion of GPF at the C-terminus of GSNOR 
driven by native promoter sequence and transformed this con- 
struct into the GSNOR null mutant hot5-2 (Lee et al, 2008). 
We carried forward three homozygous transgenic lines that 
express GSNOR-GFP (Figure 3A). Mutation of GSNOR results 
in reduced plant height and an increased number of inflo- 
rescences (Lee et al., 2008; Kwon et al., 2012). We quantified 
these differences in mature plants and found that while hot5-2 
and hot5-4 null mutants and corresponding wild-types (Col-0 
and WS, respectively) produce a similar number of first-order 
inflorescence stems (arising from the rosette), mutants produce 
two-fold more second-order branches and often produce third- 
order branches, which are not observed in wild-type plants 
(Figures 3B-D). All three complemented lines showed quanti- 
tative restoration of wild-type height and branching patterns 
(Figures 3C,D), indicating that the GFP fusion did not com- 
promise GSNOR function. We further observed reduced tri- 
chome branching in hotS null plants (Figures 3E,F), a previously 
unreported phenotype. Mutants primarily produced trichomes 
with only two branches, rather than the three or four typical 
of wild type (Figure 3F). 70% of hot5-2 trichomes were two- 
branched compared to virtually zero in Col, and similar differ- 
ences were detected between WS and hot5-4 (Figure 3G). The 
absolute number of trichomes was marginally higher in WS leaves 
compared with those of hot5-4, in agreement with Holzmeister 
et al. (2011), while trichome abundance did not differ between 
Col and hot5-2 (Figure 3E, lower panel), indicating the differ- 
ence in trichome number must be due to a modif)'ing gene 
present only in the WS background. Quantitation of trichome 
branching in the complemented lines showed that the GSNOR- 
GFP fusion protein restored branching to wild-type levels 
(Figures 3F,G). 

Based on complementation of hotS phenotypes, the GFP 
fusion seemed suitable to assess GSNOR localization. Whole 
seedlings at the cotyledon stage and flower stage 13 were 
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FIGURE 3 I Expression of GSNOR-GFP under its own promoter rescues 
multiple aspects of the hot5-2 null phenotype. (A) Immunoblots of leaf 
protein extracts probed with anti-Arabidopsis GSNOR antibodies (top) or 
anti-cytosolic GAPDH blots (middle). Rubisco large subunit (Coomassie 
stain, bottom). Plant genotypes are as indicated with CGI, 2, and 3 being 
independent homozygous 13 lines expressing GSNOR::GSNOR-GFP in the 
hot5-2 background. (B) Shoot systems of indicated genetic backgrounds 



with primary inflorescence overlaid in red. Bar: 5 cm. (CD) Shorter plant 
height (C) and inflorescence branching order (D) are rescued to wild-type 
levels in CGI, 2, and 3. (E) Trichomes with reduced branching are more 
numerous in hot5-2 and hot5-4 rosette leaves than in respective wild-type 
backgrounds. Error bars: STD (F) Left to right: trichomes from hot5-2, Col, 
Col. Bar: 0.2 mm. Red Arrows indicate trichome branches. (G) Trichome 
branching is rescued to wild-type levels in CGI, 2, and 3. 



imaged to observe the tissue and subcellular localization of 
GSNOR-GFP (Figure 4). Fluorescence was observed through- 
out seedling and floral structures (Figures 4A-C), which cor- 
relates with gene expression data in public databases (Toufighi 
et al., 2005). The seedling apical meristem and root tip exhibited 
notably intense fluorescence (Figures 4A,D), while cotyledon, 
hypocotyl, root, and petal vascular tissue signals were highest 



(Figures 4A-C). Detailed observation of sections through the 
root tip (Figures 4E-G) show distribution in all root cell types 
with diffuse cytosolic and nuclear localization, but dramatic 
exclusion from the nucleolus. This is also evident in the elongat- 
ing zone of the root, although nucleoli are considerably smaller 
in these cells (Figure 4H). GSNOR-GFP could also be detected in 
anther filaments, ovary, stigma, and petals (Figures 4I,J,M,N,P) 
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FIGURE 4 I GSNOR is expressed in various organs and developmental 
stages in Arabidopsis. Localization of GSNOR was observed in plants 
transformed with GSNOR::GSNOR-GFP in the hot5-2 background. All 
plants were homozygous for the GSNOR::GSNOR-GFP transgene, with the 
exception of those used for images in (l-K). (A) Whole seedling. (B) 
Cotyledon. (C) Flower, stage 13. (D) GSNOR-GFP distribution and 
localization in root tip cells in optical cross section through the middle of 
the root. (E) GSNOR-GFP localization in root epidermal cells. (F, G) Optical 
cross section of GSNOR-GFP localization in root cortex cells at two 
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different magnifications. Arrowheads: nucleolus. PI, Propidium Iodide 
staining; DIG, Differential Interference Contrast microscopy. (H) 
GSNOR-GFP localization in the root elongation zone. Bar: 20|j,m (D-H). 
(I-P) Localization of GSNOR-GFP in stage 14 flowers (I) stamens and 
petals (J) pollen (K) anther filaments (L) ovary (M) petals (N) seed at 
bending cotyledon embryonic stage (0) and stigma (P). Arrowheads in (K) 
denote hot5-2 pollen that do not express GSNOR-GFP due to segregation 
of the transgene in the heterozygote. Bar: 500 |im (l,J) 200 |im (IVl), or 
40 (im (K,L, N-P). 



of Stage 14 flowers and was particularly enriched in pollen and 
seed (Figures 4K,0). 

MULTIPLE PATHOGEN RESPONSE GENES ARE DOWNREGULATED IN 

hot5-2 

To assess global changes in the transcriptome due to GSNOR 
absence, we performed microarray analysis on 4-week leaves 
of hot5-2 plants grown on a 12 h light cycle since the leaves 



of Col and hot5-2 plants are most morphologically compara- 
ble at this age (Lee et al., 2008). Imposing stringent criteria 
of >2-fold changes with p < 0.01 (FDR-corrected), we found 
99 and 170 transcripts up- and downregulated, respectively, in 
hot5-2 compared to Col-0. A complete list of transcripts is pro- 
vided in Supp. Table 2 categorized by pathway as curated in 
MAPMAN (Usadel et al, 2009). This list was used to test for 
enriched categories of regulated genes (Usadel et al., 2006). 



www.frontiersin.org 



November 2013 | Volume 4 | Article 430 | 7 



Xu et al. 



GSNOR conservation and function 



One group of enriched genes was the "Stress Response" cat- 
egory (Table 1). Of 19 genes in this category, 13 were down- 
regulated "biotic stress" genes, including pathogenesis-related 
protein 1 (PRl, ~6. 9-fold), consistent with data from Feechan 
et al. (2005), four potential pathogen receptors (Atlg59218, 
At3g04210, At3g50470, At3gll010), PR5, and other predicted 
herbivore and pathogen defense proteins. Of two upregulated 
"Stress" category genes, the defensin protein PDF 1.2 was also 
reported to be upregulated in GSNOR antisense plants (Espunya 
et al, 2012). 

These results are consistent with the involvement of NO in 
pathogen responses, and prompted us to employ the Arabidopsis 
eFP browser (Toufighi et al, 2005) and TAIR gene annotations 
to determine if other pathogen response-linked genes, including 
those in the "Not Assigned" category (the largest category with 
76 genes) were also differentially regulated in hot5-2. The eFP 
data include response to the oomycete Phytophthora infestans, 
the bacterium Pseudomonas syringae (virulent and avirulent), 
the fungi Botrytis cinerea and Erysiphe orontii, as well as induc- 
tion by non-host bacteria and several elicitors. As highlighted 
in Supp. Table 2, 43 additional genes that can be linked to 
pathogen response are down-regulated (e.g., ACD6, At4gl4400; 
DIRl, At5g48485), making a total of 56 of 170 downregulated 
genes linked to pathogen response. In contrast, only 13 addi- 
tional upregulated genes can be linked to pathogen response. 
Since ~12% of Arabidopsis genes are pathogen defense-related 
(The Arabidopsis Initiative, 2000), we conclude the absence 
of GSNOR disproportionately downregulates pathogen response 
genes. 



AN UNUSUAL CLASS OF GLUTAREDOXINS ARE UPREGULATED IN 

hot5-2 

Given the significant interplay between GSNO levels and GSH- 
controUed redox homeostasis (Staab et al, 2008), it was of 
interest that genes in the "Redox" category were also enriched. 
Intriguingly, six cytosolic ROXY-class glutaredoxins (GRXs) and 
an atypical chloroplast-localized thioredoxin (TRX) were upreg- 
ulated in hot5-2 (Table 2 and Supp. Table 2). GRXs and TRXs 
are small oxidoreductases that regulate the thiol redox state of 
other proteins (Meyer et al, 2012). Four of the six up-regulated 
GRX genes are closely linked on chromosome 4 and share 91-95% 
identity, while two others on chromosomes 5 and 1 are 72-75% 
and 58-60% identical to the chromosome 4 genes, respectively, 
and 59% identical to each other. All six GRXs belong to a plant- 
specific GRX family containing a monocysteine active site. Thus, 
these GRX proteins likely act as monothiol GRXs (Herrero and 
De La Torre-Ruiz, 2007). One cytosolic, monothiol ROXY GRX 
(ROXY20), which is divergent from the upregulated genes (~40% 
identical), is downregulated in hot5-2. Examination of publi- 
cally available expression patterns of these regulated ROXY genes 
(Schmid et al, 2005) indicates that during normal plant growth 
the hot5-2 upregulated ROXY genes have low transcript levels in 
stems, senescing leaves, and the shoot apex, while ROXY20 shows 
the opposite pattern. Thus, these GRX proteins are likely to serve 
different functions. Currently the redox targets of these cytoso- 
lic GRXs and the unusual chloroplast TRX are unknown, but 
increased expression of the corresponding genes indicates that 
mutation of GSNOR alters redox homeostasis and potentially the 
targets of these oxidoreductases. 



Table 1 | Changes in "Stress Responsive" Transcripts in hot5-2 vs. wild-type. 



Gene 


Protein 


Stress 


Fold-change (log2) 


p-value 












UPREGULATED 










At3g04720 


PR4 (PATHOGENESIS-RELATED GENE 4) 


Biotic 


1.397 


0.001 


At5g44420 


PDF1.2; PDF1.2A; LCR77; Defensin 


Biotic 


1.017 


0.001 




At2gl4610 


PRl (PATOGENESIS-RELATED GENE 1) 


Biotic 


-2.78 


0 


At4g 19820 


Glycosyl hydrolase family 18 protein 


Biotic 


-2.25 


0 


At2g37570 


SLT1 (SODIUM AND LITHIUM TOLERANT 1) 


Abiotic 


-2.005 


0.001 


At1g73330 


DR4 (DROUGHT REPRESSED 4); peptidase inhibitor 


Biotic 


-1.882 


0 


At3g50480 


HR4 (HOMOLOG OF RPW8 4) 


Biotic 


-1.861 


0.001 


At1g75040 


PR5 (PATHOGENESIS-RELATED GENE 5) 


Biotic 


-1.821 


0.001 


At3g11010 


AtRLP34 (Receptor-like Protein 34) 


Biotic 


-1.8 


0.002 


At1g24020 


MLP423 (MajorLatexPprotein-LIKE PROTEIN 423) 


Abiotic 


-1.65 


0 


At3g50470 


HR3 (HOMOLOG OF RPW8 3) 


Biotic 


-1.502 


0 


At2g43510 


ATTI1; serine endopeptidase inhibitor 


Biotic 


-1.414 


0 


At3g04210 


Disease resistance protein (TIR-NBS class) 


Biotic 


-1.391 


0.001 


At2g43530 


Trypsin inhibitor Defensin-like protein 


Biotic 


-1.282 


0.001 


At1g59218 


Disease resistance protein (CC-NBS-LRR class) 


Biotic 


-1.272 


0.001 


At3g48080 


Lipase class 3 family protein 


Biotic 


-1.238 


0.001 


At1g72260 


THI2.1 (THIONIN 2.1) 


Biotic 


-1.117 


0.001 


At2g03720 


MRH6 (morphogenesis of root hair 6) 


Abiotic 


-1.083 


0.008 


At1g19670 


C0R1 (CORONATINE-INDUCED PROTEIN 1); CHL1 (CHLOROPHYLLASE 1) 


Abiotic 


-1.013 


0.001 


p < 10-^ listed as 


= 0. 
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Table 2 | Changes in "Redox" Transcripts in hot5-2 vs. wild-type. 

Gene Protein Function/Localization 



Fold-change (log2) 



p- value 



1 IDDCr^l II ATCr\ 










At4g 15690 


Glutaredoxin S5 


Glutaredoxin; cytosolic; CxxS active site 


2.788 


6.94 E-05 


At4g 15670 


Glutaredoxin S7 


Glutaredoxin: cytosolic; CxxS active site 


2.758 


8.85 E-05 


At4g 15660 


Glutaredoxin S8 


Glutaredoxin: cytosolic; CxxS active site 


2.088 


0 


At4g 15700 


Glutaredoxin S3 


Glutaredoxin; cytosolic; CxxS active site 


1.786 


0 


At1g03020 


Glutaredoxin SI 


Glutaredoxin; cytosolic; CxxS active site 


1.595 


0 


At5g04720 


ACHT5 (ATYPICAL CYS HIS 


Thioredoxin; chloroplast; CGGC active site 


1.128 


0.001 




RICH THIOREDOXIN 5) 








At5g 18600 


Glutaredoxin S2 


Glutaredoxin: cytosolic; CxxS active site 


1.064 


0.001 




At5g11930 


Glutaredoxin CIO 


Glutaredoxin: cytosolic; CxxS active site 


-1.155 


0.001 


At1g20620 


CAT3 


Catalase; peroxisome; cytosol? 


-1.108 


0.002 



p< l(r^ listed as = 0. 



CHANGES IN hot5-2 GENE EXPRESSION INDICATE OTHER PROCESSES 
IMPACTED BY NO HOMEOSTASIS 

Another enriched category of genes that presents a consistent 
picture of changes in the GSNOR mutant is the "Signaling" cate- 
gory, in which 14 of 15 genes are down regulated and half of the 
downregulated genes are involved in calcium signaling, including 
multiple calmodulins, calmodulin-like proteins and calreticulins 
(Supp. Table 2). The single upregulated gene in this category is 
ATCPl (CALCIUM BINDING PROTEIN 1, At5g49480), which 
also has sequence similarity to calmodulin. These data provide a 
direct link of GSNO to calcium signaling. 

Although "Transcription" was not a specifically enriched cate- 
gory, consideration of regulated transcription factors (Supp. Table 
1) shows that three of the most highly upregulated genes encode 
BHLH proteins involved in iron deficiency responses (BHLHIOO, 
At2g41240; BHLH039, At3g5698; BHLH038, At3g56980) (Wang 
et al., 2013). The potential targets of these transcription fac- 
tors in leaves are unknown, although in roots they can activate 
genes required for iron uptake. Among downregulated genes 
in the transcription category, it is notable that there are three 
AP2/EREBP transcription factors, which are involved in ethylene 
responses, including the most strongly downregulated gene in this 
category TINY (At5gll590; >5-fold decreased) (Sun et al, 2008). 

DISCUSSION 

The pleiotropy of Arabidopsis GSNOR loss-of-function mutants 
indicates that enzymatic control of GSNO levels is essential for 
competitive viability (Lee et al., 2008; Kwon et al, 2012) (Feechan 
et al, 2005). In vitro, GSNOR catabolizes GSNO and HMGSH, 
metabolites generated from NO- and formaldehyde, respectively. 
Because NO can lead to the generation of reactive nitrogen species 
(RNS) and formaldehyde can produce ROS (Staab et al, 2008; 
Kubienova et al, 2013), GSNOR is likely critical to regulation 
of downstream physiological and pathological effects of RNS 
and ROS in plants. The evolutionary conservation of GSNOR 
in eukaryotes and bacteria, which was recognized over a decade 
ago (Liu et al., 2001), further indicates this enzyme has similar 
functions in multiple domains of life. Here we have identified 
conserved features of GSNOR genes and their encoded proteins 



in plants, determined aspects of tissue and cellular localiza- 
tion of GSNOR, and shown that its absence impacts expression 
of pathogen response, redox, and calcium signaling genes in 
Arabidopsis. While this work is primarily descriptive, it offers 
insight into the molecular consequences of GSNOR loss and 
therefore can serve as a prospectus for future mechanistic studies. 

GSNOR was found as a single-copy gene in 35 of 41 plant 
species, with seven instances of recent duplication events (data 
not shown). De Smet et al. (2013) theorize that the preponder- 
ance of single copy genes in organisms that have undergone whole 
genome duplication (such as Arabidopsis) is not random, but 
rather indicates that these genes impair fitness when present in 
multiple copies, possibly due to overly-robust activity of their 
encoded proteins. In any event, GSNOR copy number appears to 
be under strong selective pressure. 

All putative eukaryotic GSNORs we examined show high con- 
servation and are unusually rich in cysteines, which are highly 
positionally conserved among plant orthologs (Figures 1, 2 and 
Supp. Figure 1). Common functions of protein cysteine residues 
include coordination of metal atoms (i.e., copper and zinc), cova- 
lent catalysis, extracellular adhesion, and redox sensing (Wang 
et al., 2012). Most of the ex-zinc cysteines in Arabidopsis 
GSNOR were found to be inaccessible to solvent (Figure 2; 
Crotty, 2009), suggesting their primary function may be struc- 
tural. Intriguingly, however, three cysteines that are positionally 
conserved between plants and animals are solvent accessible 
(Figure 2, blue patches). These residues may serve as sites of 
post-translational regulation via, for instance, glutathionylation 
or S-nitrosation. Residues that bind HMGSH or coordinate zinc 
(Kubienova et al., 2013) were found in all predicted plant pro- 
teins (Figure 1 and Supp. Figure 1), indicating GSNO reductase 
and HMGSH dehydrogenase activities are probably general fea- 
tures of plant GSNORs. GSNORs from Arabidopsis and tomato 
exhibit ~ 10-fold higher velocity of NADH-dependent GSNO 
reduction than NAD+ -dependent HMGSH oxidation to formyl- 
glutathione (Crotty, 2009; Kubienova et al., 2013), but the high 
ratio of NAD^" to NADH in most living cells would favor the 
dehydrogenase reaction. This could be circumvented by cofactor 
recycling, wherein NADH produced from HMGSH oxidation 
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is employed for GSNO reduction (Staab et al., 2008). Thus, 
GSNO catabolism and formaldehyde detoxification may be par- 
tially interdependent processes. This scenario is not unreasonable, 
as ROS and RNS are known to contribute to the formation of one 
another (Molassiotis and Fotopoulos, 2011). The mitochondrion 
is a major site of ROS production in eukaryotes. GSH and NO 
concentrations are reported to be high in Arabidopsis mitochon- 
dria (Wang et al, 2010; Koffler et al, 2013), and mitochondrial 
enzymes have been identified as nitrosation targets (Millar and 
Day, 1996; Palmieri et al., 2010). There is no transcriptional 
evidence for mitochondrion-localized GSNORs apart from one 
Physcomitrella homolog (Figure 1), but the existence of cryptic 
target peptide-encoding sequences in frame with GSNORs from 
diverse taxa (data not shown) underscores possible mitochondrial 
localization in a common evolutionary ancestor. 

Visualization of GSNOR-GFP provided clear evidence for 
cytosolic and nuclear localization of this protein throughout 
the plant, with potentially higher concentrations in vascular tis- 
sues and very noticeable exclusion from the nucleolus. Though 
GSNOR lacks a nuclear localization signal, it may be trans- 
ported in association with another protein. We saw no evi- 
dence of GSNOR-GFP in mitochondria, but the upstream region 
with potential to encode a mitochondrial targeting peptide was 
not included in our construct. Our results of organ and tis- 
sue localization do not fully agree with a previous report from 
Espunya et al. (2006); while they found high levels of activ- 
ity in roots, no activity was detected in the hypocotyl or 
cotyledons, which show significant fluorescence in our study. 
Consistent with our observations, their immunocytochemistry 
showed GSNOR in all cell types of the root meristematic zone. 
However, the inner cortex of the root elongation zone appeared 
to lack protein when localized by immunocytochemistry, in con- 
trast to our observations. Overall, GSNOR appears to function 
in essentially all plant cell types. These protein data are aug- 
mented by publically available transcript analysis in Arabidopsis, 
which show virtually ubiquitous expression of GSNOR mRNA, 
further supporting the role of this enzyme in multiple plant 
processes. 

A well-documented feature of GSNOR nuU Arabidopsis is a 
multi-branching phenotype (Lee et al., 2008), which has been 
suggested to arise from impact on auxin transport and func- 
tion, and/or cytokinin signaling (Kwon et al, 2012). Replacement 
of wild-type GSNOR with GSNOR C-terminally fiised to GFP 
produced plants with quantitatively normal branching patterns, 
indicating GSNOR-GFP effectively functions like the wild-type 
protein in processes required for normal branching. The GSNOR- 
GFP transgene also rescued reduced trichome branching, a pre- 
viously unrecognized phenotype of the null mutant. Arabidopsis 
trichomes are single cells, so reduced branching results from alter- 
ation of cell morphogenesis, which involves a wide range of basic 
processes (Htilskamp, 2000). How absence of GSNOR reduces 
trichome branches is not known, but notably, Arabidopsis tri- 
chomes have two to four times the concentration of GSH com- 
pared to other epidermal cells (Gutierrez-Alcala et al., 2000). 
High GSH may increase the deleterious effects of the absence 
of GSNOR specifically in these cells. GSNOR-GPF also res- 
cued this phenotype, implying a role for NO homeostasis in 



basic cell morphological processes. The pleiotropy of the GSNOR 
loss-of-function phenotype suggests enhanced nitrosation of 
one or more proteins interferes with their normal physiolog- 
ical functions. However, there are only a few known in vivo 
targets of nitrosation in plants: the Arabidopsis cytokinin sig- 
nal relay kinase AHPl is negatively regulated via nitrosation 
of a conserved cysteine (Feng et al., 2013), while nitrosation 
of the F box protein TIRl has been proposed to modulate 
auxin signaling (Terrile et al., 2012). Thus, the strong mor- 
phological phenotypes of hot5-2 and hot5-4 plants may reflect 
uncoordinated cytokinin and auxin crosstalk caused by aber- 
rant nitrosation of TIRl and AHPl. Observed changes in 
pathogen sensitivity may have a similar etiology. Indeed, sali- 
cylic acid (SA) signaling has been shown to be downreguated 
by GSNO-mediated nitrosation of NPRl (Tada et al, 2008). 
Further work with the GSNOR null mutant will be required 
to determine any alterations in nitrosated or glutathionylated 
targets. 

Comparison of transcript levels between a GSNOR null 
mutant and wild type Arabidopsis on microarrays revealed 
an enrichment in transcripts involved in pathogen responses, 
redox regulation, and calcium signaling, as well as upregula- 
tion of transcription factors involved in iron responses and 
downregulation of ethylene responsive transcription factors. 
The extensive downregulation of pathogen defense related and 
responsive genes supports the observations of Feechan et al. 
(2005), who reported the nuU mutant had reduced R gene medi- 
ated, basal and non-host resistance to pathogens. While the 
Arabidopsis homolog of human CRl was not transcriptionally 
higher in the hot5-2 mutant, as might be expected if it pro- 
vided compensatory activity, GRXs and a TRX were notably 
upregulated. TRXs and GRXs can mediate denitrosation and 
deglutathionylation, respectively (Benhar et al, 2009; Zaffagnini 
et al., 2013), such that they could function to reverse an increase 
in these modifications due to excess GSNO in the mutant. 
Although the contribution of TRXs and GRXs to the acclima- 
tion of GSNOR-deficient plants to nitrosative stress has not been 
assessed, these microarray data strongly implicate the involve- 
ment of GRXs in plant SNO homeostasis. The microarray data 
also support proposed linkages between NO and calcium sig- 
naling (Courtois et al., 2008), as five calmodulins or calcium 
binding proteins are downregulated and one is upregulated in 
hot5-2. 

We also compared the results of our microarrays with previ- 
ous studies aimed at identifying NO regulated genes that used 
NO donors or NO synthase inhibitors to induce changes in 
gene expression in Arabidopsis (Parani et al., 2004; Besson-Bard 
et al., 2009). We found very limited overlap with potentially 
NO-regulated transcripts. For genes induced by application of 
0.1 or 1.0 mM SNP (Parani et al., 2004), only four genes 
were also upregulated in hot5-2: the calmodulin-like ATCPl 
(At5g49480), WRKY40 (Atlg80840), a UDP glucosyl trans- 
ferase (Atlg05560) and a 67 amino acid unknown (At4g27654). 
Compared to transcripts regulated by a NOS inhibitor, in 
which downregulated genes were proposed to be normal tar- 
gets of NO, again only four transcripts behaved similarly 
in our experiments (DINIO, At5g20250; XTR8,At3g44990; 
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LTPl,At2g45780; unknown, At3g56360). We may see little overlap 
with these prior investigations because they examined short-term 
manipulation of NO levels, whereas GSNO and NO are chroni- 
cally deregulated in hot5-2. In any event, changes in gene expres- 
sion in the absence of GSNOR reflect a long-term metabolic 
adjustment required to cope with excess GSNO and the pathways 
it normally regulates. 

GSNOR clearly plays a role in biotic stress adaptation, but 
how a GSNOR-DOD fusion such as that identified in a homolog 
from apple (Supp. Figure 1) would function in vivo is not 
clear, and this may be an artifact of early phase sequencing. 
Using the Toronto Bio-Analytic Resource (BAR) Arabidopsis 
gene expression data compendium (Toufighi et al., 2005), it 
was revealed that GSNOR transcription strongly correlated with 
NINJA [r^ = 0.62, a jasmonic acid (JA) response corepressor] 
and PMR5 (r^ = 0.59, a protein whose absence affords greater 
resistance to bio trophic fungi that cause powdery mildew). 
Both PMR5 and NINJA negatively regulate pathogen-induced 
defense signaling — PMR5 contributes to JA-independent fungal 
disease susceptibility (Vogel et al., 2004), while NINJA activity 
is curtailed following JA-induced, COIl -dependent proteasomal 
degradation of JAZ repressor proteins (Pauwels et al., 2010; 
Sheard et al, 2010). While the association of these pathogen 
response genes with GSNOR is only correlative, it can be 
inferred that GSNOR might also work to dampen biotic stress 
responses in the absence of elicitation. This would explain why 
GSNOR over-expression and RNAi-mediated knockdown served 
to respectively diminish and enhance systemic acquired tolerance 
and basal tolerance to a _R syringae and Peronospora parasitica 
(Rusterucci et al., 2007). This also harmonizes well with the 
observation that GSNOR positively affects SA signaling (Feechan 
et al., 2005), since JA and SA operate antagonistically to one 
another. 

In summary, GSNOR appears to be a ubiquitously-expressed, 
cytosol-localized protein that regulates shoot morphology, 
pathogen defense responses, and NO homeostasis. Aberrant 
nitrosation of auxin, cytokinin, and SA response regulator pro- 
teins, among others, likely contribute to aspects of hotS null 
mutant phenotypes. Diminished branching in GSH-rich tri- 
chome cells further underscores the role of GSNOR in main- 
taining the cellular reduction potential, and its conserved, 
solvent-accessible cysteines may function as NO sinks or serve a 
regulatory role. The upregulation of transcripts of a class of GRXs 
is a particularly promising discovery, as some GRXs and TRXs 
catalyze deglutathionylation and denitrosation. Understanding 
how GRXs may compensate for loss of GSNOR and how GSNOR 
activity may be regulated through its conserved, solvent accessible 
cysteine residues will help to clarify the role of GSNOR in plant 
biology. 
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